Astronomy & Astrophysics manuscript no. bergvall' 1303 


©ESO 2013 


March 8, 2013 





Lyman continuum leaking gaiaxies - 
Searcli strategies and iocai candidates'^ 

Nils Bergvall^, Elisabet Leitet^, Erik Zackrisson^, and Thomas Marquart^ 

^ Department of Physics and Astronomy, Uppsala University, Box 515, SE-751 20 Uppsala, SWEDEN 
e-mail: nils .bergvall@physics .uu. se , elisabet . leitet@physics.uu.se , thomas.marquart@physics.uu.se 

^ Department of Astronomy, Stockholm University, Oscar Klein Center, AlbaNova, Stockholm SE-106 91, Sweden 
e-mail: ez@astro .su.se 

Submitted 

ABSTRACT 

Context. Star-forming dwarf galaxies leaking Lyman-continuum (LyC) radiation may have played an important role in the reionization 
of the Universe. Local galaxies exhibiting LyC leakage could shed light on the escape mechanisms, but so far only two such cases 
have been identified. Here, we argue that this meager harvest may be caused by unsuitable selection criteria, one of which is a high 
equivalent width in Ha emission. Such a criterion will guarantee a high production of LyC photons but also indicates a high column 
density in the neutral gas, eff'ectively ruling out LyC escape. 

Aims. We want to investigate whether the lack of local LyC emitters can be caused in part by biased selection criteria, and we present 
a novel method of selecting targets with high escape fractions. By applying these criteria, we assemble a sample of observation targets 
to study their basic properties. 

Methods. We introduce a new selection strategy where the potential LyC leakers are selected by their blue colours and weak emission 
lines. We also take a closer look at the properties of 8 LyC leaking candidates at z ~ 0.03, selected from the Sloan Digital Sky Survey 
(SDSS) which we observe with ESO/NTT in broadband B and Hof. 

Results. We find that 7 of the 8 target galaxies are involved in interaction with neighbours or show signs of mergers. In 7 cases the 
young stellar population is clearly displaced relative to the main body of these galaxies. In about half of our targets the absorption 
spectra show young post-starburst signatures. The scale lengths in Ha are typically 30% smaller than those of the stellar continua, 
which is characteristic of galaxies influenced by ram pressure stripping. We tentatively identify a few conditions favourable for 
leakage: 1) the combined eff'ects of ram pressure stripping with supernova winds from young stars formed in the front, 2) merger 
events that increase the star formation rate and displace stars from gas, 3) starbursts in the centres of post-starburst galaxies, and 4) a 
low dust content. 

Key words. Galaxies: intergalactic medium - starburst - fundamental parameters - evolution. Cosmology: diff'use radiation. 
Ultraviolet: galaxies 



1. Introduction 

After cosmic recombination the "dark ages" followed. It is an 
epoch during which no astronomical light sources existed and 
the gas in the universe remained neutral. A few hundred mil- 
lion years later, the universe underwent reionization. The seven- 
year WMAP data suggests that this process started at z > 
10 (Komatsu et al. 2011), and observations of Gunn-Peterson 
troughs in the spectra of high-reds hift quasars places the com- 
pletion of this process at z ~ 6 (e.g. Fan et al.|2001 2006 ). The 
sources responsible for reionization still, however, need to be 
identified. 

Accreting black holes and population III stars could have 
contributed to the LyC budget during the reionization epoch (e.g. 



Madau et al. 2004| 


Trenti & Stiavelli|2009l|Volonteri & Gnedin| 


2009( [Mirabel et a] 


J2011), but galaxies of low to intermediate 



masses seem to be the top candidates. Star formation models 
and high-z observations (e.g [Stark et aL]|2007[ |Bouwens et al. 



201 1| ) imply that galaxies were already in place at z ^ 10 and 
that starburst dwarf galaxies may have played a dominant role. 
Cosmological simulations give further support to the importance 



* Based on observations made with ESO Telescopes at the La Silla 
or Paranal Observatories under programme ID 080.B -0899(A) 



of galaxies during the reionization epoch, even though predic- 
tions vary as to the exact ma ss scales of the objects giving the 
larges t LyC contribution (e.g|Gnedin et al.||200"8{ [Wise & Cen| 



[20091 [Razoumov & Sommer-Larsen|2010[|Yajima et al.|2011 

Most attempts to detect the escape of hydrogen-ionizing (i.e. 
Lyman continuum, hereafter LyC) radiation from local galaxies 
have failed (e.g. |Leitherer et al.|p^95l [Deharveng et al.||200H 
Grimes et al.|2009|). Even though a few galaxies, in cluding our 



own Milky Way (B land-Hawthorn & Maloney|1999|), display in- 
direct signatures of leakage ( [Heckman et al.|2Ql HlZastrow et al. 
201 1 ), the starburst galaxies Haro 11 (Bergvall et al. 2006; Leitet 
,et al.,201 i nand Tol 1247-232 (Leitet et al. 2012, A&A, submit- 
ted) are the only cases from which LyC photons have been di- 
rectly observed. The question remains, however, to what extent a 
complete census of LyC leakage in the low-redshift Universe has 
yet been reached. Local LyC leakage candidates are typically se- 
lected on the basis of having strong emission lines of high excita- 
tion. While this guarantees that the target galaxies are producing 
ionizing photons in large numbers, this strong-line criterion will 
also bias the selection towards galaxies with a high H i content 



The 



Bergvall et al. | |2006| ) detection was later disputed by 



Grimes 



|et al.| ( [2007| ), but the recent re-analysis by [Leitet et al. | ( 20TT] ) confirms 
the detection and places the leakage at a level of /esc ~ 3% 
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relative to the ionizing flux. In this case the H ii regions will be 
radiation bounded, which actually prohibits LyC escape. Here, 
we explore an alternative strategy for selecting potential LyC 
leakers from the Sloan Digital Sky Surve}j^ (hereafter SDSS) 
in the hope of increasing the number of low-redshift galaxies 
with known LyC escape. Detailed studies of such local, spatially 
resolved targets may then shed ligh t on the mechanisms respon- 
sible for LyC escape (e.g. Haye s et al.|20Q7 ) in a way that would 
be staggeringly difficult for the LyC leakers identified at higher 
redshifts. Eventually, this should have important implications for 
our understanding of cosmic reionzation. 

A flat ACDM cosmology with Ho=70 km s"^ Mpc"\ 
11^=0.3, and flo=0.1 was adopted throughout the article. All 
magnitudes are given in the AB system, and flux densities in 



2, The escape fraction 

Direct measurements of leaking LyC radiation from the reioniza- 
tion epoch galaxies are not feasible, since the Lya forest is too 
dense at these redshifts. Instead we are limited to the study of 
analo gous objects at z < 4 - 5 ( |Kim et al.|2QQ7[|lnoue & Iwata 
|2008| ). The amount of ionizing radiation leaking from galaxies 
can be quantified by the LyC escape fraction. For local galaxies 
the absolute escape fraction (/esc) is mostly used, defined as the 
ratio of the observed flux density of LyC photons (fLyCohs) over 
the intrinsic LyC flux density produced by the stars (fLyC,mt)- The 
parameter fiyCM is usually derived from the extinction corrected 
Hof flux. 

For distant galaxies the Ha line is shifted out of the optical 
wavebands, and fiyCM is therefore difficult to access. Instead the 
relative escape fraction (/esc,rei) is often used, defined as the ratio 
of the intrinsic Lyman break amplitude derived from SED mod- 
els, (fi5oo/fLyc)mu ovcr the observed amplitude, (/i5oo//L3;c)obs, 
corrected for the IGM optical depth. To compare the relative es- 
cape fraction of distant galaxies with the absolute escape fraction 
of local galaxies, the former needs to be corrected for the internal 



1.00 



absorption by dust at 1500 A, A 1500: 

^sc — ^sc,rel X 10 = 
(/l500 //LyC )int _ . IGMx ^ 1 a-0.4-Ai5oo 



(1) 



(/i500//l3/c) 



obs 



X expir'l^y^') X 10" 



where rj^^^ is the IGM optical depth for the LyC photon. 

fiyCohs is usually measured at rest wavelength ^^900 A, but 
for z ~ 1 galaxies the LyC has be en sampled at shorter wave- 
lengths even down to ^700 A (Malkan et al.||2003 



2007) [2QT0t [Cowie et al.|[2009t [Bridge et al.]|201Q| ). In the lit- 



Siana et al. 



erature, the intrinsic amplitude has often been assumed to be 
(/i50o//L}/c)int= 1-2.5 (3-7 in units of fy) for a normal stellar pop 
ulation based on Starburst99 (e.g. Siana et al.|2007 



et aL 



Iwata et al. 



2009 ). This agrees with predictions from our models ( [Zackrisson 
]2QQT] ) 



It has been suggested that th e ionizing emissivity of galax- 
ies is evolving with redshift (e.g. |Inoue et a l. 2006 ; Siana et al. 
12010'). Many investigations at z 3 do indicate high escape 
fractions, both for single detections and stacked samples, while 
the low-z universe shows a different picture. In the following we 
have made an effort to summarize and visualize the existing data 
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Fig. 1. The absolute escape fraction, /esc, as a function of red- 
shift. The data points are a mixture of single detections with 
errors, detections in stacked samples with errors and upper lim- 
its from stacked samples. The absolute escape fraction has been 
calculated from reported (fisoo / fiyc) ohs values, using a homog- 
enized treatment of the IGM opacity and the intrinsic ratio 
(fi50o/fLyc)mu as described in the text. 



in one single plot. To do this we have obtained the observed ra- 
tios (fisoo / fiyc) ohs from the literature. Since the treatment of the 
IGM opacity varies between differe nt authors, we have homog- 



http ://w w w. sds s . org/ 



enized r^^^ using the prescription in Nestor et al.| (|201 1 ) around 
z ~ 3, and Siana et al. ( 2010 ) for z ~ 1. I his means, tor example, 
that a galaxy at z=3.2 will have tJ^^^=0.95 while one at z=1.3 

will have t^^^=0.59. One should be aware that the dispersion 
in o ptical dept h can b e quite substantial for different sightlines 
(Inou e & Iwa ta 2008 ; Nestor et al.|2011 )), but here we can only 
apply the mean. We have also used a single value for the in- 
trinsic Lyman break amplitude (fi5oo/fLyc)mu even though this 
is known to be strongly dependent on parameters such as the 
burst age and metallicity, and could in principle also be evolv- 
ing with redshift. Here, we have used the intrinsic ratio derived 
from Haro 11, which is the only LyC leaking galaxy yet where 
the SED model has been constrained by observations in multiple 
filters ranging from the far-ultraviolet ( FUV) to the near inf rared 
(NIR), including the Ha emission line ( Hayes et al.|2007 ). The 
intrinsic ratio for Haro 11 was found to be (/i5oo//L3;c)int=l-5 
(Leitet et al. 2012, A&A, submitted), somewhere in the middle 
of the range of values used in previous papers. The value 1.5 
is valid for a stellar population older than ~20 Myr, using the 
Zackrisson et al. ( |2001| spectral synthesis model. For the galax- 
ies where the LyC is sampled at shorter waveleng ths, we have 
made a correction to the intrinsic ratio following |Siana et aL 
('2007'). The internal absorption, E(B-V), was adopted as either 
the value reported by the author, or where no value was reported, 
as the mean value from similar galaxies at the same redshift. To 
derive the extinction at 1500 A, we have used the ,Calzetti| ( | 1 99l\ 
extinction law where Ai5oo=10.33xE(B-V). 

The escape fraction as a function of redshift is shown in 
Fig. [T] The dispersion in the (fi5oo/fLyc)mt model data has not 
been included in the error budget, only the observational uncer- 
tainties. The data are obtained from the following papers. At 
z ~ 3 we use the stacked sample from |Steidel et al.| p001| ), 
the stack ed sam ple and the one confirmed single d etection from 
Shap ley et al.| ([2006), the stacked sample from Iwata et al. 
f2009| ), the stacked sample and the z=3.795 single detection 
(with a lower limit dependent on the IGM opacity) from Vanzella 



2 



Bergvall et al.: Lyman continuum leaking galaxies 



et al. ( 2010a ), the stacked samples of Lyman break gal axies 
(LBGs) and Lya emitter s (LAEs) fromj Nestor et al. (20111, and 
the stacked sample fro m|Boutsia et al.[ (2011). At z ^ 1 we us e 
the stacked sample s of|Malkan et al.|(|20 03 ); Sia na et aL |(|2007|); 



Cowie et al. (2009 ); Siana et al. 



2010| ) and Bridge et al. 



(2010). 



In the local universe we use the single detection from [Leitet 
|et al.| ( |20TT] ). In the local universe we have also included a new 
detection, for the galaxy Tol 1247-232, as well as the stacked 
spectrum of 9 FUSE galaxies (Leitet et al. 2012, A&A, submit- 
ted). We have chosen not to include upper limits based on single 
galaxies. 

From the limited information we have, the results suggest an 
overall lower /esc for the present day galaxies than for the galax- 
ies at z 3, but the relation is weakened with the recent results 
from |Vanzella et al.| (|2010a') and |Boutsia et al.|p01 It where only 
upper limits could be derived for their stacked samples. There is 
also a substantial risk that some of the z ^ 3 detections are actu- 
ally low redshift interlopers ( V anzella et al.|2010b|[20T2| ). 

Here, we argue that the few reports of LyC leaking galaxies 
at lower redshifts, and the apparent decrease in /esc at z < 3, 
may not only be the result of an evolving escape fraction. There 
are also several possible biases to consider. One such bias con- 
cerns at what wavelengths the Lyman continuum is sampled. All 
investigations at z ~1 have been carried out in the rest wave- 
length range 700-850 A, while most investigations at z ~3 and 
in the local universe sample around 900 A. Even a small amount 
of dust, allowing some LyC photons to escape at 900 A, would 
have an increasing effect with decreasing wavelength. The gen- 
eral importance of dust and how it may contribute to the de- 
tection rates over a large redshift range is further discussed in 
Sect. |5.2[ Another possibility is that if there would still exist pop 
III stars at z ~1, the nebular "Lyman bump" emissio n just be- 
low the Lyman limit ( Inoue||2010[ Inoue et aL]|2011| ) would be 
missed. 

Another important factor is the bias caused by selection ef- 
fects. The main difference between surveys at z ~3 and those 
performed in the local universe is the target selection. Most sin- 
gle detections have been made at z~3 when narrowband imaging 
is utilized just below the Lyman limit in multi-object fields. The 
selection criteria for these objects are only weakly constrained 
by including LBGs and LAEs within a certain field and redshift 
range. On the contrary, the selection of galaxies observed in the 
local universe is heavily constrained by the assumed properties 
of LyC leakers. 

One urgent question now emerges - Are we really probing 
the local equivalent to the LyC-emitting galaxies found at high 
z? Here technology is the limiting factor: multi-object fields can- 
not be observed, and we can only use spectroscopy (currently 
HST/COS). The only thing we can improve is the selection strat- 
egy of the targets. In the following sections we will present a 
novel method to select galaxies that are the most likely to have 
high escape fractions of ionizing photons, and take a closer look 
at a selection of z ~ 0.03 LyC leaking candidates. 

3. How to optimize the search for leal<ing galaxies 

3.1. Basic conditions for LyC leakage 

In the local universe, dwarf galaxies involved in intense star- 
bursts have been the main targets in searches for leaking galax- 
ies. These galaxies have strong emission-line spectra with Ha 
equivalent widths typically in the range EW(Hck) = 500 - 1000 
A. In a galaxy with a significant leakage however, the nebular 
contribution would be much lower (oc(l-fg^c) in the picket-fence 



model, see below), due to the lack of available gas to ionize. 
Although the colours would be blue, such a galaxy would have 
a 'dull' optical spectral appearence with weak emission lines. 
Therefore, they could naively be regarded as uninteresting and 
escape selection. If the leakage is due to a global ionization of 
a significant fraction of the Hi envelope (in the idealized case 
a truncated Stromgren sphere) the degree of ionization should 
be high. If, on the other hand the radiation is leaking through a 
cleared channel, the so-called 'picket-fence' model, it should be 
more normal. From what we know, these aspects of the expected 
observational properties of a LyC leaking galaxy have been ne- 
glected in previous searches. In this paper we describe a method 
to select galaxies with these pr operties. Our method finds further 
support from recent results in Nestor et al. ( 2011| ) of a weaker 
hya emission from LyC leaking galaxies than non-leakers. 

An additional selection criterion to be added is the age. The 
early stages of the reionization probably takes place in the ap- 
proximate redshift range 9<z<ll. This means that the starburst 
must be younger than the corresponding timespan, i.e. <100 Myr. 
We thus have to limit our models discussed below, such that the 
maximum age of the young component is < 100 Myr. 

3.2. Dust attenuation 

An important issue concerns how to deal with the presence of 
dust in leaking galaxies. In the previous section we mentioned 
two basic ideas related to the conditions for LyC escape - the 
density bounded model (truncated Stromgren sphere) and the 
picket-fence model. Between these two extremes we find the 
more realistic density bounded ionization cone model where part 
of the light from the galaxy is escaping through the cone and part 
is absorbed in other directions. In the picket-fence model it is as- 
sumed that the radiation is leaking out through gas- and dustfree 
holes in the interstellar medium while the non-leaking regions 
are dusty. This condition may develop when gas outflows from 
the star forming regions open channels towards the halo. To cre- 
ate a density bounded cone model, the Lyman continuum flux 
needs to be intense enough to ionize all atomic hydrogen gas 
along the cone, allowing excess photons to disperse into the in- 
tergalactic medium. In this case the dust should be evenly dis- 
tributed over the source, in the simplest case like a foreground 
screen. In the haloes of the galaxies the radiation density is much 
lower than in the starburst region. Under such conditions, even 
though the gas is ionized, small dust particles may absorb LyC 
photons without being destroyed (Boulanger et alT]|1988| ). The 
critical radiation density is ~ a few eV cm"^ . The Hor fluxes we 
derive from the target galaxies are typically 2-3xlO^^W, corre- 
sponding to a LyC flux of 2-3xl0^^eV s"^ Assuming the pro- 
duced photons to spread out isotropically, we obtain a radiation 
density of leV cm"^ at a distance of ~ 3 kpc from the source. 
Since this is below the critical radiation density for grain de- 
struction, we would expect dust to survive at larger distances and 
contribute to attenuation of the light. A simila r condition seems 
to be at hand in e.g. M82 ( Telesco et al.|1 989 ). The actual situ- 
ation in leaking galaxies is probably somewhere in between the 
picket-fence and truncated Stromgren sphere models and in the 
following we will consider both of them to fully constrain the 
models. 



3.3. A tentative z -0.03 sample 

In order to have a closer look at galaxies with the desired proper- 
ties in the local volume, we performed a first selection of leakage 
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Table 1. Basic data for the SDSS sample. Columns: 1) Running number 2) Name of the target galaxy 3) Right ascension 4) 
Declination 5) Redshift 6) Distance, obtained from NED 7) Absolute g magnitude based on SDSS model magnitudes 8) Model u-g 
colour 9) Model g-r colour 10) Fiber u-g 11) Fiber g-r 



No 


Name 


RA (a) 


Dec (6) 


z 


r 


M, 


u-g 


g-r 


(u-g)/ 


(g-r)/ 






(J2000) 


(J2000) 




Mpc 












(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


1 


J010724.76+133209.2 


01'^07'"25' 


+ 13^32'"09^ 


0.038 


154 


-19.3 


0.72 


0.20 


0.69 


0.15 


2 


J024352.54-003708.4 


02M3^53^ 


-00^37'"08' 


0.030 


120 


-17.5 


0.32 


-0.15 


0.56 


-0.01 


3 


J025325.28-001356.6 


02^53^25^ 


-00^13'"57^ 


0.027 


109 


-18.5 


0.67 


-0.17 


0.33 


-0.09 


4 


J075313.34+123749.1 




+ 12^37'M9^ 


0.029 


123 


-17.6 


0.74 


0.12 


0.69 


0.11 


5 


J080754.64+141043.8 


08^07'"55^ 


+ 14^10'«44^ 


0.029 


122 


-19.3 


0.81 


0.29 


0.59 


0.13 


6 


J085642.06+123158.3 


08^56'"42^ 


+ 12^3r"58^ 


0.029 


127 


-18.3 


0.82 


0.24 


0.66 


0.21 


7 


J100712.21+065735.6 




+06^57'"36^ 


0.032 


136 


-18.9 


0.61 


0.19 


0.70 


0.05 


8 


J101007.57+033130.5 


10H0^08^ 


+03^3r"3r 


0.032 


137 


-18.3 


0.65 


0.25 


0.71 


0.17 



Table 2. Derived properties. Columns: 1) Running number 2) D4000n index 3) E(B-V) derived from the Ha/Hfi ratio after correction 
for underlying absorption 4) Fiber u-g colour, corrected for extinction assuming a foreground screen 5) Fiber g-r colour, corrected 
for extinction 6) The [Oni]/15007/HyS ratio 7) Oxygen abundance derived from the ON and ONS methods according to Pilyugin] 
et al.| ( [2QTQ] ) 8) The exponential scale length in the Ha continuum 9) The Ha scale length 10) Equivalent width of Ha in emission 



derived from our measurements of the SDSS spectra (NEW) and as given in the SDSS data bank (SDSS) 11) The Ha equivalent 
width within a circular aperture of 10 arcseconds radius. 



No 


D4000n 


E(B-V) 


(U-g)/,c 


(g-r)/,c 


[Oiii]/H^ 


12+log(0/H) 


^Ha-cont 


hna 


EW(Ha) (A) 


EW(HQf) (A) 












NEW/SDSS 


ON/ONS 


kpc 


kpc 


NEW/SDSS 


Total 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


1 


1.01 


0.15 


0.44 


-0.02 


1.5/1.9 


8.42/8.50 


1.91 


1.41 


67/67 


81 


2 


0.91 


0.04 


0.50 


-0.05 


2.0/2.2 


8.34/8.46 






166/180 


66 


3 


0.82 


0.16 


0.11 


-0.28 


3.8/4.0 


10.2/8.7 


1.56 


1.14 


230/228 


170 


4 


0.95 


0.22 


0.32 


-0.15 


3.1/3.4 


8.11/8.28 


1.36 


0.59 


146/146 


78 


5 


0.95 


0.13 


0.37 


-0.03 


1.8/2.0 


8.38/8.45 


1.30 


1.04 


119/118 


70 


6 


0.98 


0.11 


0.47 


0.08 


2.1/2.6 


8.11/8.25 


2.32 


1.32 


79/79 


72 


7 


0.94 


0.11 


0.53 


-0.08 


3.0/3.3 


8.14/8.29 


1.44 


1.07 


100/108 


72 


8 


1.02 


0.15 


0.45 


-0.02 


2.0/2.4 


8.06/8.40 


1.61 


1.60 


71/70 


51 



candidates in the SDSS at z -0.03 using the following selection 
criteria: 

- 0.025 <z< 0.040 

- 20 A< EW(Ha) < 200 A 

- 2.5 < Ha/H/3< 4 

- EW(H/?) > 10 A 

- 05/H/3 > 2 

- u-g < 0.7 (prior to extinction corrections) 

In addition to this, AGNs were deselected according to the 
SDSS classification and the criterion that non-AGNs should have 
a Gaussian width of Ha of (T//Q,<5x(z-hl) A. Any remaining 
AGNs were later efficiently removed using the BPT diagnos- 
tics ("Baldwi n et al.||1981| ). The lower redshift limit was chosen 
to make sure that the galaxies were distant enough for a large 
fraction of the light from the target to enter the 3" aperture of 
the SDSS spectrograph. The limiting magnitude for the SDSS 
spectroscopic survey is m^ -17.5. This corresponds to ~ -18 
at z=0.03. To account for a reduced amount of nebular emis- 
sion if the galaxy is leaking, we selected galaxies in the range 
20A<EW(HQf)<200A. The lower limit was chosen in order to 
avoid too large uncertainties in the extinction corrections. The 
upper limit in combination with the colour criterion, u-g<0.7, 
strongly increases the probability of finding leaking galaxies in 
the sample, although it does not completely exclude non-leaking 



cases. We will argue below that a minimum requirement to ful- 
fil the properties of three of our five reddest target galaxies is 
at least fesc =20%. This is valid under the assumption that the 
galaxies are dust- free. But as we will show below, if we take the 
dust attenuation into account, we arrive at a minimum absolute 
escape rate of -5%. Galaxies with lower escape rates will be dif- 
ficult to identify due to the uncertainties in the model predictions 
and we therefore chose this value, corresponding to EW(Ha)~ 
200A, as an upper limit. Our Ha observations described below, 
show that this criterion holds globally for all of the eight galax- 
ies finally selected for observations. The extinction corrections 
are based on the Ha/H/3 ratio and at low EW:s the reliability of 
the correction is very sensitive to how one applies the correction 
for underlying absorption. To avoid galaxies with high extinction 
we set an upper limit to Ha/Hj3 of 4 and to avoid problems with 
correction for underlying absorption we demand EW(H/3) to be 
larger than lOA. The 05/HJ3 criterion guarantees that the exci- 
tation is high and thus that the light is dominated by massive 
shortlived stars. According to our models, this criterion is valid 
about 20 Myr after an abrupt shutdown of a starburst. The colour 
criterion is based on model predictions of a young stellar popu- 
lation mixed with an old host and includes a normal amount of 
interstellar extinction. The galaxies were also selected to be visi- 
ble from the ESO telescopes. Of the -350000 objects catalogued 
within this redshift interval in the full SDSS data release 7, only 
230 galaxies met the criteria (<0.1%, see Sect. 5.2 for a discus- 
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between the predicted values for this galaxy. Typical abundances 
are 12+log(0/H) ~ 8.3, corresponding to 40-50% solar metal- 
licity (usin g the solar oxygen abundance as derived by |Asplund| 



4000 4500 5000 5500 6000 6500 
Wavelength (A) 



Fig. 2. Optical spectra of the target galaxies, obtained from the 
3" diameter fiber of the SDSS. The spectra have been shifted ver- 
tically with an arbitrary amount in order to make the individual 
spectra more clearly visible. The numbers refer to the running 
numbers in Table [T] 

sion). Next we calculated the extinction corrections based on the 
Ha/ll/S ratio, after correction for underlying absorption, using a 



method based on the D4000n break developed by |Mattsson & 



Bergvall (2009). The eight most promising LyC leaking candi- 



dates were then singled out for more detailed investigations. The 
SDSS spectra are displayed in Fig.|2] 

The basic properties of the target galaxies are listed in 
Table fTl and the derived properties in Table [2] In column 7 of 
Table [2] we give two alternative determinations of the oxygen 
abundances. Accurate abundances can be calculated if we know 
the electron temperature of the ionized gas. Normally the weak 
[O in]/14363 line is used for this purpose. If the line is too noisy, 
semiempirical methods based on strong emission lines have fre- 
quently bee n used, the most popula r ones being the R23 (|Pagel| 
|eraLlp^79l ) and the N2 methods ( [Denicolo et"ari[2QQ2| . The 
reliability of the methods have to be tested from time to time 
as more accurate data become available. The strong line meth- 
ods are therefore continuously being improved. A discussion of 
these issues with regard to SDSS data is found in the paper by 
[Shi et al.| ( |2Q10| ). Here we have used an extension of these meth- 
ods as described by P ilyugin et al.| ( |2010} . These are based on a 
combination of oxygen and nitrogen lines (ON method) and oxy- 
gen, nitrogen and sulphur lines (ONS method). One should take 
these abundances with a pinch of salt since, if the galaxies are 
leaking, the H 11 regions are not radiation bounded and the equa- 
tions used to derive the oxygen abundances are not fully appli- 
cable. Anyway, the derived abundances agree very well between 
the two methods and also between the eight galaxies except for 
galaxy no. 3. This galaxy is attributed a higher metallicity than 
the others, outside the metallicity range used for calibration of 
the ON and ONS methods. This probably explains the deviations 



et al.|2004| ). This (40% solar) is the metallicity we have adopted 
for our model comparisons in Fig.|4](see below). Comparing the 
eight potential LyC leakers with other local star forming galax- 
ies, the galaxies of our sample are residing among the most lu- 
minous BCGs in a metallicity-luminosity plot (see Fig. |3]). But 
one should remember that the completeness limit at the typical 
redshift is ~ -18 so potential leakers may exist also at lower 
luminosities. The metallicities of the targets used in the diagram 
are mean metallicities as presented in Table |2] Data for the other 
galaxies are collected from various sources (Bergv all et al., in 
preparation). The known weak LyC leaker Haro 1 1 ( [Leitet et aL 



[201 1 ), often used as a local analogue of LBGs at z ~3, is located 
to the extreme right in the distribution (see Fig.[3|. 

The low internal reddening estimated from the Ha/H/S ratio 
indicates that the weak Ha emission is not due to any abnor- 
mally low dust content for these galaxies. Fig. [4] shows how a 
random sample of galaxies at z -0.03 with 20 A<EW(HQf)<200 
A is distributed in the u-g/g-r diagram after correction for dust 
attenuation. The mean attenuation for the comparison sample in 
Fig.|4]is A(V)~0.45. Some guidance to how the colours relate to 
the spectral features are obtained from Fig.|5] We note that Ha is 
still in the r window, i.e. not redshifted out. If the emission lines 
and the nebular continuum are weak, the colours essentially de- 
pend on the stellar continuum and this is the major trend we see 
in the evolutionary track with the highest escape fraction. The 
dramatic change in colour for this track is due to the change in 
age of the stellar population with the youngest population having 
the bluest colours. 

In Fig. [4] we also show the predicted colours of star form- 
ing galaxies fulfilling the EW(HQf) criterion, but with diff'erent 
amount of leakage, /^-parameters and age. The models consist of 
mix a young and an old component, including a nebular compo- 
nent, and are described in more detail in Appendix A. In the fig- 
ure we also see the positions of our selected target galaxies, be- 
fore and after corrections for dust attenuation, assuming a fore- 
ground screen. Unfortunately, the two confirmed leakers Haro 1 1 
and Tol 1247- 2321.4 are not observed in the Sloan survey so we 
cannot compare to their positions in the diagram. 

The total escape fraction for galaxies fulfilling our selec- 
tion criteria can be estimated if we know the true dust attenu- 
ation correction in the UV. This correction is difficult to obtain 
so we will discuss two extreme cases - the picket-fence model 
and the homogeneous dusty foreground screen model. The ac- 
tual conditions in the galaxies are probably intermediate be- 
tween these cases. We split the escape fraction into two compo- 
nents: fesc=fesc,gas'fesc4usu whcrc fesc,gas IS the fraction escaping 
through the neutral hydrogen gas and fescdust is the fraction re- 
maining after passing through the dust. Let us first have a look 
at the picket-fence model. Here it is assumed that the leaking ra- 
diation is unattenuated so fesc,dust=^- To allow for a comparison 
between observations and the models in Fig. [4] we should how- 
ever correct the non-leaking regions for dust attenuation, using 
the Ha/U/S ratio. We do not know the exact proportions of the 
dust attenuated part and the leaking part. However, the model 
can help us to make a fair guess. Most of the targets have colours 
that, before correction, closely agree with the model shown with 
a dotted line. This model has a leakage of fesc= 0.40. The so- 
lution is not unique. In fact, acceptable models can be found in 
the range 0.20< fesc <0.60. But the 40% model gives the best 
agreement with the observed colours and EW(HQf) of the tar- 
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Fig. 3. The luminosity-metallicity diagram for the z -0.03 
galaxies. In the plot other local star forming galaxies of types 
Irregular, Blue Compact (BCG) and Low Surface Brightness 
(LSBG) are shown (Bergvall et al. 2012, in prep.). The eight 
potential LyC leakers, marked with large squares, are distributed 
among the most luminous BCGs. Also shown are two known 
local LyC leakers, Haro 11 and Tol 1247-232. We have also in- 
dicated the much studied low metallicity galaxy I Zw 18. Typical 
Icr errors are indicated to the upper left. 



gets. If we accept a deviation between model and observations of 
30% in EW(Ha) and 0.05 magnitudes in the colours, this model 
satisfies 80% of the galaxies in the upper envelope. So we will 
assume that the escape rate is at least 40%. If we add dust to 
the 60% non-leaking part, the total colours will become redder. 
Therefore the leakage must be more than 40%. By gradually in- 
creasing the fraction of leakage and correcting the non-leaking 
part for dust attenuation we find that the colours (the squares in 
the diagram) will experience a small shift towards bluer colours 
with the reddest colours corresponding to an absolute leakage of 

fesc~fesc,gas ~50±10%. 

Returning to Fig. |4j we look at the second option - a dusty 
foreground screen. To correct the optical data for dust attenua- 
tion we used the observed Ha/H/? ratios compared to the nor- 
mal recombination ratio to derive the extinction coefficients. In 
Fig.|4]the corrected, dust- free colours are shown as triangles. The 
positions of the triangles can now be directly compared to the 
predictions from the dust-free model. We could of course make 
individual model fits to each of the target galaxies and then de- 
rive the predicted escape fractions from the observed EW(HQf). 
However, because the of the simple assumptions we make about 
the dust attenuation, w hich may vary with po sition and age of the 
stellar population (e.g. Chariot & Fall 200 ()|), we prefer to look at 
trends. Obviously most of the data points agree with a model that 
has leakage of /esc = /escgas ~ 0.95 and /esc,dust=l. Now we have 
to find the true value of /escdust for a foreground screen with the 
derived mean extinction coefficient. Using the relations derived 
by ,CalzettietaL| ( ,1994 ) and Cartle dge et al. ( 2005 ) we find that 
< /escdust >=0.1. Thus the absolute escape fraction as measured 



in front of the screen is /esc ~ 0.95 -0.1 ~ 0.1. A more conserva- 
tive value of the minimum escape rate of the sample galaxies is 
estimated to be /esc ~0.05. Thus it is clear from this discussion 
that our target galaxies have properties that according to simple 
models allow for a substantial (5-50%) leakage of LyC photons. 



Fig. 4. Comparison between model predictions and observed 
colours. The observed target colours are marked with squares. 
The triangles show the colours after correction for dust attenua- 
tion assuming a homogeneous foreground screen. The open cir- 
cles show a random sample of 1000 star forming SDSS galaxies 
at z=0.03 with 20<EW(HQf)<200 A, also corrected for dust red- 
dening. All corrections are based on the Ha /HyS ratio and are rel- 
evant for redshift z=0.03 and 40% solar metallicity. Typical er- 
rors in the colours are 0.03-0.04 magnitudes but some of the out- 
liers may be faint galaxies with large errors. Overlaid are tracks 
showing Zackrisson et al. (2001) model colours of galaxies with 
and without LyC leakage. Only parts of the evolutionary tracks 
where the EW criterion is fulfilled are shown. We look for cases 
which can be approximated by a mixture of a starburst with con- 
stant star formation rate and an old component. The duration of 
the starburst is 10 Myr. These tracks can be characterized by 
two parameters: 1) the birthrate parameter /?, i.e. the ratio be- 
tween present star formation rate and the mean past star forma- 
tion rate, 2) /esc- The age of the old component was assumed to 
be 10 Gyr. The lines represent the following scenarios: Hatched: 
b=l, /esc=0 ; Dotted: b=2, /esc=0.4; Solid: b=30, /esc=0.95. The 
temporal evolution is from the top of the diagram. In the upper 
tracks we notice a kink on the blue side. This marks the end of 
the starburst. The nearly horizontal part of the two upper tracks 
correspond to the post- starburst phase. 



3.4. Optical Imaging 

3.4.1 . Observations and data reduction 

The z -0.03 galaxies were imaged with the EMMIR and SUSI2 
at the ESO New Technology Telescope (NTT) at La Silla, Chile, 
December 13-15, 2007. The three nights were dark and had good 
photometric conditions. The average seeing was 1.2''- 1.5'', with 
occasional periods of ~ 2" seeing. Ha images were obtained with 
EMMIR. We used the narrowband filters #599 and #600 (both 
with FWHM ~ 71 A), having central wavelengths correspond- 
ing to the redshifted Ha line of the galaxy. The Hof offline fil- 
ter (#768) with central wavelength 6014 A was selected to only 
include the stellar continuum not too far from the emission line. 
Images in the Cousins B band were obtained with SUSI2. The in- 
strument set-up and integration times can be seen in Table[3] The 
data were reduced with the NOAO/IRAF package, using stan- 
dard procedures including bias subtraction, flat-field correction, 
spectrophotometric standard stars, and sky subtraction. The Ha 
images were produced from the online images after subtracting 
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Fig. 5. The spectrum of target galaxy No. 1 with the SDSS ugr 
filters overlaid. 



the offline images, scaled to make the stars in the field disappear 
after subtraction. A comparison with the SDSS EW(HQf) was 
carried out as a final check. The total EW(Ha) within a radius of 
10 arcseconds is given in Table 2. These values are all within 
our search criteria. In Table [3] are also given the exponential 
scale lengths in Ha and the Ha-continuum. As can be seen from 
Fig. |7] the Ha-emission is close to spherical, or the morphol- 
ogy is so distorted that it is difficult to define a symmetry plane. 
When we derived the scale lengths we therefore did not correct 
for inclination eff'ects. We integrated radially and azimuthally 
from a position corresponding to the brightest part of the Ha im- 
age. In the derivation of the scale length, the central regions were 
excluded and only the exponential part of the images outside the 
very centre was used. 



4. Results 

4.1. Morphologies 

The target galaxies are shown in Fig.[6j [7] and [8] The first figure 
shows the isodensity contours of the galaxies in the Ha contin- 
uum band. The second figure shows the Ha emission in isophotal 
contours with Ha continuum in greyscale. The third figure shows 
the EW(Ha) greyscale map overlaid on the isodensity contours 
of Fig. [6] From a look at Fig. [6] we can see that it appears that 
there are two types of galaxies. One is dominated by a nearly 
spherical component (no. 1 & 8). Others show a disturbed mor- 
phology wit h a a few galaxies (no. 4 , 5 & 7) displaying a tadpole 
structure ( va n den Bergh et al.|1996| ) in Ha. Although our galax- 
ies are fairly bright, albeit sub-L*, this is uncommon among nor- 
mal galaxies at these luminosities. The general impression is that 
most of the target galaxies appear to be dynamically disturbed, 
either with clear signs of interaction or clearly non-symmetric 
light distribution. Two galaxies, no. 1 and 8, at first appear to be 
in equilibrium, but also in these cases the Ha maps overlaid with 
the continuum maps reveal clear asymmetric structures. 



4.2. Young stars and ionized gas ■ 
leakage 



implications for LyC 



How can we use our observational data to better understand the 
conditions for leakage? We require from the global data that 



EW(HQf) should be low while at the same time the stellar popu- 
lation should be young. But how can we find a stronger support 
for leakage locally? We show two examples of our targets in 
Fig [9] The figure shows maps of a colour index and EW(Ha). 
The index is defined by the ratio of the ffux in the B band and 
the Ha continuum ffux as measured in the EMMI fflter #768 at 
the eff'ective rest wavelength 6026 A, which at the redshift of 
our galaxies corresponds to a rest wavelength of 5850 A. We 
call this index B-EMMI768. It is a strong indicator of young star 
forming regions. In both examples shown in the ffgure, we see 
that the young population is forming a 'channel' from the cen- 
tre out to the very limit of the optical image. The conditions for 
leakage will be higher the closer that young stars are to the halo 
region. If the gas is ionized by stars at the edge of the main body 
(deffned by the noise limit of the optical image), the observed 
EW(Ha) should be low. This is what we observe in galaxy no. 7 
on the left. We can derive a rough estimate of the mean density 
in the ionized gas from 



eff 



471 



,0 n 

J -ex 



mi,R])e 



-lr{l,R)lhHu 



dl (2) 



where R is the projected distance from the centre (sometimes 
called the impact parameter), I is the distance in the tangential 
direction along the line of sight and r = ^/(P' -\- R^) is the ra- 
dial distance. We assume a spherical distribution of the ionized 
gas and an exponential density decrease where hnii is the scale 
length in mean density of the ionized gas. j^Ha is the surface 
brightness in Ha. a^^. is the eff'ective recombination coefl&cient 
ofHjS (3.034 lO'^^crrPs'^) and r^^ is the em ission line ratio be- 



tween Ha and HJ3 (2.85), both obtained from Hu mmer & Storey 
(1987 ) assuming an electron temperature of Te=lO'^K and a den- 
sity of 100 cm"^. Both these parameters are only weakly depen- 
dent on densities at low density and can be applied here even if 
we ffnd much lower densities, h is the Planck constant, vh^ is 
the frequency of the HJ3 emission line, rip^o and rie^ are the cen- 
tral proton and electron number densities of the ionized gas in 
cm"-^ (riefi ~ ripfi [= hq] is assumed here), f(r) is the density de- 
pendent fflling factor, i.e. the fraction of space that is fflled with 
ionized gas, assuming a bimodal density distribution. 
From this we obtain 



no = 



1/2 



(3) 



and n = noe~^^^"" at / = 0. To proceed, we need to know 
the scale length of the ionized gas and how the fflling factor 
varies with distance from the centre. Here we will assume that 
the mean gas density follows the luminosity distribution in the 
Ha continuum. From Table [2] we see that the scale length in 
Ha is shorter than in the Hof-continuum which is as expected 
since the Ha emissivity goes as n^. The calculations are not very 
sensitive to the scale length we choose but more to the change of 
fflling factor with density. Determinations of the fflling factor in 
the Milky Way and the Andromeda galaxy ( B erkhuij s en||2004[ 
|Berk huijsen e t al.,2006 ) show a nearly linear dependence of the 
fflling factor with the reciprocal electron density, / oc n~^. The 
fflling factor thus increases with decreasing density and in the 
diff'use ionized gas (DIG) it is fairly high. In regions of active 
star formation, it may be more poro us due to the winds form the 
star forming regions ( Martin I99l\ |de Avillez & Breitschwerdt 
2004 ). Here we will assume that the fflling factor follows this be- 
haviour and increases with distance from centre. For simplicity 
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Fig. 6. Images in the Ha continuum of the target galaxies obtained with EMMI/NTT at ESO. East is left and north is up. The contour 
levels differ by a factor of 2. The numbers refer to the running numbers in Table [T] 
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Fig. 7. Images of the target galaxies. East is left and north is up. The contours show the Ha emission, while the inverted grey scale 
background shows the continuum from the Ha off-line spectral region. The numbers refer to the running numbers in Table [T] 



we will assume that the filling factor is constant along the line 
of sight towards the faint blue regions since the distance from 
the centre, and thus the filling factor, is almost the same along 
a scale length distance from the densest region in the tangential 
direction. We will have a look at two extreme values of the fill- 
ing factor, f=0.0 2 and f=0.2, defining a range typically found 
in the DIG (e.g. Miller & VeiHeuxl[jn03| [Martel et al.|[2Q04l 
[Berkhuijsen & Muller||2008 ). We will assume an electron tem- 



perature of 10"^ K and that the scale length of the mean density 
of ionized gas is the same as for the Ha continuum, hnii =1.5 
kpc. Using /^Ha = 3.8-10"^erg s"^ cm"^ sr"^as measured towards 
the optical edge of the blue region of galaxy No. 7, we obtain a 
mean density of the ionized gas at about 4 kpc fom the centre of 
< He >~ 0.01-0.04 cm"^. This is about 2 dex lower than the den- 
sities in the largest H II regions in other galaxies, having sizes of 
100-200 pc or more (e.g. Scoville et al.|[2001| ) and a few times 
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Fig. 8. Here is shown the Ha equivalent width in grey scale with the Ha continuum background shown in the isophotal map. The 
flux levels are arbitrary and are intended to give a good overview of the stellar + free-free emission component. The numbers refer 
to the running numbers in Table [T] 

Table 3. The NTT z -0.03 sample observational details. Column: 1) Running number 2) Name of the target galaxy 3) Ha online in- 
strument/filter (exposure time in seconds) 4) Ha offline instrument/filter (exposure time in seconds) 5) the B images instrument/filter 
(exposure time in seconds). 



No 


Object 


Ha online 


Hot offline 


B filter 


(1) 


(2) 


(3) 


(4) 


(5) 



J010724. 
J024352. 
J025325. 
J075313. 
J080754. 
J085642. 
J100712. 
J101007. 



76+133209.2 
54-003708.4 
28-001356.6 
34+123749.1 
64+141043.8 
06+123158.3 
21+065735.6 
57+033130.5 



EMMIR/#600 (3100) 
EMMIR/#599 (3100) 
EMMIR/#599 (3100) 
EMMIR/#599 (2190) 
EMMIR/#599 (2190) 
EMMIR/#599 (2100) 
EMMIR/#599 (4380) 
EMMIR/#599 (2190) 



EMMIR/#768 (1770) 
EMMIR/#768 (1770) 
EMMIR/#768 (1770) 
EMMIR/#768 (930) 
EMMIR/#768 (930) 
EMMIR/#768 (900) 
EMMIR/#768 (1860) 
EMMIR/#768 (930) 



SUSI2/#811 
SUSI2/#811 
SUSI2/#811 
SUSI2/#811 
SUSI2/#811 
SUSI2/#811 
SUSI2/#811 
SUSI2/#811 



(800) 

(1200) 

(1200) 

(960) 

(960) 

(960) 

(960) 

(960) 



lower than the DIG in the halo regions in normal disk galaxies 
(e.g. [Collins & R^[2QQH [Miller & VeiTIeuxl[2QQ3| . The low 
density we find is more similar to that of the outflow regions 
of starburst galaxies, e.g. M82 ( [Bauer et al.|2Q07[ ). Although the 
uncertainties in our derived numbers may be significant, the very 
low values lends support to favourable conditions for leakage. 

We can also imagine another scenario where the ionizing ra- 
diation is not produced in the outskirts. Maybe there is a more 
centrally located ionization source that is so powerful that it can 
ionize the medium all the way into the halo in a certain direc- 
tion. If the scale length of the ionized gas is larger than that of 
the stars, EW(Ha) can (but need not) increase as we approach 
the halo and we can observe extremely high values at faint lev- 
els. We see this happen in the second example, galaxy no. 3, 
to the right in Fig [9] In this transition region we find extremely 
high values of EW(Ha), about 1000 A. This could indicate that 
the gas is ionized far outside the region where the ionizing ra- 
diation is formed, in the low density halo. Within the ionized 



medium there may be low density channels through which the 
radiation easily can escape. Now we will show below that the 
scale length in Ha in fact is systematically smaller than the stel- 
lar scale length so how can we explain the observations? If the 
galaxy is subject to ram pressure stripping, the stellar disk may 
be truncated in the direction of the movement through the inter- 
galactic medium ( [Koopmann et al.[[2006| ). We will discuss ram 
pressure stripping below. If this is the case in galaxy no. 3, the 
rapid increase in EW(HQf) is understood. The most important 
conclusion here however is that the ionizing stars in the two 
cases seem to form a 'bridge' between the main body and the 
halo that would favour leakage. 

Not all galaxies in our sample show a strong off'-center star 
formation. This could partly be due to projection eff'ects but 
also because we may be dealing with diff'erent types of galax- 
ies, where the leakage becomes possible for diff'erent reasons. 
In two cases (galaxies 1 & 8) the galaxies have very regular 
structures, deviating for the other non-equilibrium cases (see 



9 



Bergvall et al.: Lyman continuum leaking galaxies 




Fig. 9. The target galaxies no. 7, J100712.21+065735.6 (left pair) and no. 3, J025325.28-001356.6 (right pair). The left image in 
each pair shows the colour distribution based on the ratio between the fluxes in the B and Ha -continua, index B-EMMI768. Blue 
regions map young stars and red regions map old stars. The image on the right in each pair shows a map of EW(HQf), ranging 
between and 200 A in the galaxy on the left and between and 600 A in the galaxy to the right (white regions are higher values). 
Blue regions have low equivalent widths, red have high equivalent widths. The images have been selectively (depending on the 
local S/N) filtered with a median boxcar filter, using the MIDAS task filter/smooth. North is up, east is to the left. The fieldsizes are 
24"x26" and 36"x48" respectively. 



Fig|6]). Interestingly, these two galaxies are similar to the "dom- 
inant central objects" of the Lyman-break ana logues defined 
as LyC leakage candidates by |Heckman et al. ( 201l| . One of 
these galaxies. No. 1, shows a young post-starburst signature, 
i.e. strong Balmer absorption lines. As is shown in Fig.[TO| stel- 
lar post-starburst spectra, mixed with emission line spectra, are 
seen in 3 cases. The post-starburst signatures start to develop a 
few 10 Myr after the onset of a starburst and, if the starburst de- 
cays, achieve a maximum after a few 100 Myr and are gone after 
not much more than 1 Gyr. The conditions for leakage would be 
improved in such cases if the starburst that happened some time 
ago was strong enough to expel and heat the gas in the main body 
and open channels towards the IGM through which the ionizing 
photon s could escape. With v elocities of the order of a few 100 
km^-"^ ( Heckman et al.|l990| the winds need ~10 Myr to clear 
the central few kpc:s. During this epoch the 'post-starburst' pop- 
ulation is just about to emerge while hot stars are still present 
and can leak out ionizing radiation. We expect to detect the post- 
starburst population at a younger age than normal if the emission 
lines are weakened due to leakage. 

In Fig. [8] we display EW(HQf) maps. What is characteristic 
of all but one galaxy (no. 8) is the chaotic structures in EW(HQf), 
strongly deviating from the distribution of the red stellar contin- 
uum and in most cases also the Ha emission. It essentially re- 
flects the irregular distribution of the stellar population of difl'er- 
ent ages. Apparently, most of the galaxies are not in equilibrium. 
The deviation from equilibrium favours leakage since young star 
clusters may be able to separate from the gas and during brief 
period be located in regions where the column density of the ab- 
sorbing neutral gas in some directions may be significantly lower 
than normal. This is likely to occur if the dynamical instabilities 
cause the clusters to migrate outwards on short timescales or if 
the galaxy is exerted to ram pressure stripping. 

4.3. Ram pressure stripping 

A few of our targets are residing in small groups of galaxies in 
obvious interaction with the neighbours, and one may consider 
the importance of ram pressure stripping as the galaxy moves 
through the common hot intergalactic medium. Models of strip- 



ping indicate that star formation can be increased significantly 
during the process ([ Kronberger et al.|[20b 8 ; Mast ropietro et al.| 
[2009). On a fairly short time scale, -100 Myr, t he gas on the 
side facing the wind may be severely depleted ( Kronberger et al. 
|2008 ) while strongly ionizing stars are still alive in the low den- 
sity wake immediately behind. It thus seems possible that both 
the gas removal and the star formation triggering by the motion 
through the intergalactic medium could favour leakage. Does it 
have any significant eflfect in groups? Certainly the Ha morphol- 
ogy of some of our targets remind of the cometary or tadpole 
( [van den Bergh et aL 1996 ) structure of wakes predicted from 
models of g alaxies being stripped (see e.g. Fig 11. in Kapferer| 
let al.|2009| ). 

In a recent paper, Freeland & Wilcots pOl l) find that groups 
of galaxies of the size about a magnitude more massive than 
the local group can contain intergalactic clouds with densities 
10"^ < n <10~^ and velocity dispersions of a few 100 kms"^ In 
these cases, the ram pressure is sufficiently strong to completely 
remove the gas from galaxies in the mass range ~ 10^"^ At© but 
is probably insufficient for masses above 10^ At©. As discussed 
below, our galaxies are in the intermediate range and may be 
severely affected. Two cases where ram pressure effects seem be 
active, as seen from their morphology, were discussed above and 
displayed in Fig.jp] In an interesting comparative study of Virgo 
vs. field galaxies, Koopmann et al. ( 2006) find that the ratio be- 
tween the Ha scale lengths and the R-band scale lengths of spiral 
galaxies were systematically diff'erent between the two samples. 
Field galaxies had a mean ratio of /z//q,//z/?= 1.14+0.07, while the 
Virgo galaxies had an Ha scale length which was 20% lower 
than the R-band scale length. They interpret these results as an 
environmental effect. In the Virgo cluster, ram pressure stripping 
is efficiently removing the gas from the galaxies as they move 
through the cluster. Some of the Virgo galaxies show a strong 
displacement of ionized gas in a direction opposite to the centre 
direction. These also show an increased star formation activity 
on the side facing the centre. A similar process is going on in the 
Large Magellanic Cloud. 

Can we apply the results obtained by Koopmann et al. on 
our sample? The sample differs in mean luminosity but about 
half of their galaxies have absolute magnitudes in the range - 
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Fig. 10. Spectra of the three galaxies showing post-starburst 
spectra. The numbers to the right refer to the running identifi- 
cation numbers given in the Tables. The strong B aimer lines in 
absorption are indicated. 



18.5>Mfi >-19.5, which encompasses about 50% of our sample. 
Koopmann et al. also concluded that they found no trends in 
hHa/hR with neither luminosity nor Hubble type. Most of our 
galaxies show disk structure but due to the fact that they tend 
to appear disturbed we refrain from making a proper bulge-disk 
separation and we assume that we view the systems face on. We 
choose the brightest part in the Ha continuum as the centre and 
then derive the luminosity profiles in circular strips. Then we in- 
spect the profiles and fit a linear relation in R-ji to the nearly ex- 
ponential part. In all cases we leave out the central region, corre- 
sponding to a radius of 3-4''. In Table 2 we list the scale lengths 
in Ha and the Ha continuum, excluding no. 2, which is of an 
M51 companion type. The continuum wavelength (~ 6014A) is 
fairly close to the R band. The R band also contains Ha but in 
our case the contribution is less than 10% to the R-band so we 
can more or less directly compare our scale lengths to the results 
by Koopmann et al.. In all cases the Ha scale length is shorter 
than the continuum scale length. In the mean, < hHa/hna-cont > 
= 0.71, i.e. the Ha scale length is 30% smaller. Thus the eff'ect is 
stronger than in the Virgo galaxies and signals that ram pressure 
stripping and possibly star formation triggering may be active. 
We find this to be an additional support for leakage in our galax- 
ies. 



5. Discussion 

5.1. Local versus global leakage 

In current simulations of Lyman continuum leakage from galax- 
(e.g. 



les 



Gnedin et al. 2008 [ Razoumov & Sommer-Larsen 



2010| ), the escape fraction is typically defined at the virial radius 



Tvir of the dark matter halo. Since rvir is much larger than the 
visible extent of any galaxy, this leads to a slight problem when 
attempting to indirectly assess the escape fraction from obser- 
vational data such as ours. There is simply no way that we can 
trace the light profiles all the way out to ryir, and if the broadband 
colours and Ha equivalent widths were very diff'erent once inte- 
grated out to this radius, this would indicate that the leakage is 
local rather than global. In other words, Lyman continuum pho- 



tons may be escaping from the central regions of our targets, yet 
be absorbed before reaching the intergalactic medium. 

To estimate the virial radii of our objects, we assume a stellar 
mass-to-light ratio in the range Atstars/^^ ~ 0.01-0.1 (in solar 
units) for the sample (a s der ived by mixing burst and post-burst 
populations as in Sect. 3.1), and a gas mass fraction of fg^^ ~ 
0.5. After allowing for a small amount of extinction (A(g) < 1.0 
mag), this implies a mass in detected baryons of ~ IO ^'^Mq for 
our targets. Using the baryonic Tully-Fisher relation ( |McGaugh| 
|et al.||2010 ^, we can then estimate the total dark matter masses 
of our objects. We arrive at AI500 ~ 10^^"^^ At©, where M500 
represents the mass inside which the average density is 500 times 
the critical density of the Universe. This can be converted into a 



virial mass and a virial radius using re lations in |Bryan & Norman 
(1998) and^K omatsu & Seljak| ( [2QgT] ). The result is that the virial 
radii of our target galaxies are likely to be ryir ^ 50-100 kpc at 
z ~ 0.03. In terms of angular sizes, the virial diameters are likely 
to subtend angles of ^ 3-6' in the sky at z ~ 0.03. 

The SDSS apertures (3" across) used when selecting our 
Lyman continuum leakage candidates are admittedly much 
smaller than this, and although our 5-band and Ha images can 
trace the luminosity profiles somewhat further out, the signal is 
inevitably lost in the sky noise at yUg ~ 28.6 (3cr arcsec"^), cor- 
responding to typical radii of not more than a few tenth of an 
arcminute. However, the shape of our surface brightness profiles 
(hna ~ Ikpc) suggest that the flux from the faint outskirts that 
we miss contributes no more than 10% to the total. Hence, the 
LyC escape fractions inferred by our method likely represent the 
global values, and should therefore be directly comparable with 
the results from simulations if there is negligible absorption in 
the halo. 



5.2. Evolution with redshift 

As we pointed out in Sect. 3.3, only ~ 0.1% of the galaxies at 
the redshift we are exploring fulfil our selection criteria. If these 
criteria constitute a necessary condition for a significant leakage 
it means that leakage from such galaxies is completely unim- 
portant for the ionization balance in the local universe. If the 
same conditions hold at high redshifts the leaking Lyman pho- 
tons would contribute insignificantly to the reionization. How 
do we expect the conditions to change with redshift? As the 
morphology of our targets tell us, there is a tendency for dis- 
turbed morphologies with indications of mergers or interactions 
with neighbouring galaxies or the intergalactic (intragroup or in- 
tracluster) medium. Some show tadpole structures, morpholo- 
gies interpreted as a result of mergers or ram pressure sweep- 
ing. How frequent are these objects? Locally, -0.2% of the Kiso 
UV excess gala xies (Takase & Miyauchi-Isobe|1993) show tad- 
pole structures ( [Elmegreen et al.|2012| ), i.e. the same amount 
as the galaxies that obey our search criterion. The number of 
tadpole galaxies increases dramatically with redshift and i n the 
UDF the frequency among bright galaxies is ~ 10% (Elmegreen] 
|et al.|2012j . The same authors find that the increase in SFR in the 
tadpoles is even more dramatic, of the order of two magnitudes 
from z=0 to 3, for the mass range we are dealing with. This is 
much higher than the increase in relative gas mass fraction. An 
increase in gas mass fraction clearly makes it less favourable for 
LyC escape but this may be compensated for by the increase in 
SFR in the sense that the fraction of ionized gas and the im- 
portance of SN superwinds increase. If the lifetime of the event 
when the leakage can take place is short (e.g. < a few 10 Myr 
during a stripping process), then the detection rate of the cases 
where the leakage is high increases with the volume density of 
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the objects leaking. So the increased escape fraction can partly 
be explained by the observed increase of galaxies exposed to 
ram-pressure stripping and mergers and the increased SN activ- 
ity. From these considerations we would expect a gradual change 
in LyC escape with redshift. In Fig. [T] we see a dramatic increase 
in the LyC escape fraction with redshift but at z ~ 1 we only have 
upper limits. Why is it so? One possible explanation, related to 
technical aspects, was discussed in Sect. |2] We also noted the 
possibile importance of dust attenuation. Although the dust at- 
tenuation of the galaxies at redshift of 0.03 is not 'serious' in 
the optical region, the LyC may be attenuated with up to a fac- 
tor of 10. Studies of the evolution of the dust content in the uni- 



verse g ives interesting results. Hay es et al.|p011| ) and Garel et al. 
2012| ) have investigated the correlation between Lya escape and 



dust content with redshift and find a tight relation showing that 
Lya escape increases as the dust content in galaxies decreases 
with redshift. This is of particular interest since Lya and LyC es- 
cape sh ould be related if the l eakage is from a density bounded 
region. Cucciati et al.| ( |2012| ) find that the dust attenuation in- 
creases from zero redshift, reaching a peak at z ~1 and then 
decays to a level similar to today at z ~3. All these eff'ects com- 
bined can help to explain the observations summarized in Fig.[T] 
The importance of sub-L* galaxies in the reionization process is 
further strengthened by the fact the dust content is a steep func- 
tion of stellar mass, preferentially allowing more generously for 
leakage in low mass galaxies ( ,Sawicki|2012| ). 

6. Summary and Conclusions 

In this paper we suggest that the lack of detection of leaking 
local galaxies could, at least partly, be due to a selection bias. 
Leaking galaxies, by definition, should use less photons for ion- 
ization than normal star forming galaxies and therefore the emis- 
sion lines should be weaker. However, all searches so far have 
been focused on galaxies with strong emission. Here we select 
galaxies which have fairly weak B aimer emission lines but ex- 
tremely blue colours, typical of a galaxy dominated by young 
stars but where the ionizing photons are leaking out. Using our 
spectral evolutionary model we select galaxies from the SDSS 
with the properties of a leaker. We select eight galaxies at a red- 
shift of z = 0.03 for imaging in broadband B and Ha. Although 
the sample is small we find some characteristic properties that 
deviate from those of normal star forming galaxies and support 
LyC leakage. These are 

- Among the eight target galaxies, six show clearly disturbed 
morphologies, signs of interactions, mergers or ram-pressure 
stripping. Close neighbours that also show morphological 
distorsions are found in the environments of four- five galax- 
ies. 

- In most of the galaxies we find an off'-center located popula- 
tion of young stars that, in a colour map, seems to connect to 
the halo. This shows that the conditions for leakage are im- 
proved since the hot stars are located where the gas column 
density is low. In all cases the Ha scale length is significantly 
(30%) smaller than the scale length of the stellar continuum. 
This is an indication of ram pressure stripping of the gas and 
explains our estimated low density in the ionized gas and 
some morphological features. 

- Young post-starburst signatures are seen in the stellar spectra 
of three of the eight galaxies. This may partly be an artefact 
caused by a strong leakage that makes the spectrum look as 
if the starburst has ceased. It indicates however that a major 
star formation epoch started more than ~ 10 Myr ago. During 



this epoch winds from massive stars have had time to clear 
the sight in the central region. 

In summary, we find that a plausible explanation for why 
there are so few detections of Lyman continuum leaking galax- 
ies in the local universe is that we have been ignorant of their 
presence in our search. If we instead of looking for leakage from 
strong starburst galaxies more focus on galaxies with blue con- 
tinuae and weak emission lines, our investigation of this small 
template sample shows that we may be on the right track. In a 
truly leaking galaxy, the radiation manages to reach outside the 
virial radius. Although we cannot map the whole volume out to 
this radius, our estimate is that the amount of radiation absorbed 
outside what we can observe is insignificant and that the galax- 
ies most likely are leaking. There seems to be specific conditions 
at hand that favour leakage. These are the involvement in merg- 
ers and interactions with neighbours, interaction with the inter- 
galactic medium via ram pressure stripping and a recent major 
starburst as a precursor to the present young leaking stellar pop- 
ulation. 

An additional factor that seems to play a major role over cos- 
mic time scales is the dust content. Lya emission has been found 
to anti-correlate with an increasing dust content in galaxies over 
time. The redshift dependence of LyC leakage indicates that dust 
may also contribute significantly to the possibilities for the LyC 
photons to escape. 
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Appendix A: Spectral modelling 



The m odels are based on spectra from the [Zackrisson et al. 
( 2QQ1| ) spectral evolutionary models library. The stellar popu- 
lations have a Salpeter initial mass function and a mass range of 
0.08-120 Mq. The metallicities range fr om 5% solar to solar. A 
nebular component based on CLOUDY ( |Ferland|1996[ Ferland 
et al.|1998] ) is also included, making it possible to model galaxies 
with different amounts of leakage. In the basic setup the models 



comprise a set of spectra of stellar populations of different ages 
and modes of star formation histories. The mode can be either 
constant star formation over a certain time period or an expo- 
nentially declining star formation S FR oc where t is the 
age and r is the decay timescale. 

In the modelling of the SDSS spectra in the present investi- 
gation we assume that the stellar population consists of a mix- 
ture of a young and an old population as if a burst of star for- 
mation recently took place in an old passive galaxy. To mimic 
the conditions at reionization, the young population cannot be 
too old so we assume a burst duration of 10 Myr and then an 
abrupt shut down of star formation. The metallicities are 40% so- 
lar. We mix these components in different proportions and only 
consider the case of constant star formation for the young com- 
ponent. The spectrum can be described as as a sum of young 
and old stars plus a nebular component, where the young stel- 
lar component and the nebular components are age dependent: 
^^(0 = FA,young(t) + F^^oid + F^A,Hiiit). The old component is as- 
sumed to have experienced a short initial burst of star formation 
10 Gyr ago, followed by passive evolution. We chose this age 
to agree with the cosmic look back time at peak cosmic SFR 
(Hopkins & Beacom 2006 ). Observations of local dwarf galax- 
ies also support the formation of most of the stars at z>l ( Weisz] 
|et al.||20lT] ). The star formation activity at a specific age may 
conveniently be expressed by the b - parameter, the ratio of the 
present SFR and the mean past SFR, preceding the start of the 
burst: b = SFR/<SFR>. 

The young and old components are mixed in different pro- 
portions and the relative mass fraction of young and old stars as 
function of age of the young component can be expressed as 



-^^young J ^ young 



^old 



told 



(A.1) 



where tyoung and toid are the ages of the young and old popu- 
lation respectively. We calculate a grid of models where we vary 
two parameters - the /^-parameter and the escape fraction, f^sc- 
In the code, escape fractions above fesc = can be simulated 
by numerically reducing the contribution from the nebular com- 
ponent with a suitable amount. Emission lines and the nebular 
continuum will be affected the same way. If Fnao r^fes to the 
Ha flux when there is no leakage of ionizing radiation and Fna 
is what we measure when there is leakage we find 



fesc - (1 - Fna/FHao) - (1 " ^ A,nebC I ^ A,nebCQ) 



(A.2) 



where Fx,nebc and F^^nebCo is the flux density of the nebular 
free-free and free-bound continuum emission with and without 
leakage. Consequently we may derive the equivalent width of 
Ha from 



EW(Ha) 



^Haoi^ - fesc) 



A,old 



+ F. 



A,young 



+ FA,nebCoi\ " fesc) 



(A.3) 



where F^^oid and F^^young are the flux densities of the con- 
tinuum under under Ha with contributions from the old stars 
and the young stars respectively. When we follow the evolution 
of the spectrum with time, the galaxy will experience a period 
where 20A<EW(HQf) < 200A. The larger the escape fraction, 
and assuming the same /^-parameter, the lower will be the age 
of the young population during the epoch that this criterion will 
be fulfilled. When one considers the characteristics of spectra of 
leaking galaxies one should remember that, since Fx,nebc domi- 
nates the continuum under Yia in young starbursts and Ha scales 
linearily with F^^^nebc^ EW(HQf) will initially drop only slowly 
with increasing fesc- 
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